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ABSTRACT 

We report on our early photometric and spectroscopic observations of the extremely luminous Type 
II supernova (SN) 2008es. With an observed peak optical magnitude of my — 17.8 and at a redshift 
z = 0.213, SN 2008es had a peak absolute magnitude of My = —22.3, making it the second most 
luminous SN ever observed. The photometric evolution of SN 2008es exhibits a fast decline rate 
(~0.042 mag d"^), similar to the extremely luminous Type ILL SN 2005ap. We show that SN 2008es 
spectroscopically resembles the luminous Type ILL SN 1979C. Although the spectra of SN 2008es 
lack the narrow and intermediate-width line emission typically associated with the interaction of a 
SN with the circumstellar medium of its progenitor star, we argue that the extreme luminosity of 
SN 2008es is powered via strong interaction with a dense, optically thick circumstellar medium. The 
integrated bolometric luminosity of SN 2008es yields a total radiated energy at ultraviolet and optical 
wavelengths of ^10^^ ergs. Finally, we examine the apparently anomalous rate at which the Texas 
Supernova Search has discovered rare kinds of supernovae, including the five most luminous supernovae 
observed to date, and find that their results are consistent with those of other modern SN searches. 
Subject headings: supernovae: general — supernovae: individual (SN 2008es) 



1. INTRODUCTION 

Wide-field synoptic optical imaging surveys are con- 
tinuing to probe the parameter space of time-variable 
phenomena with increasing depth and temporal cover- 
age (e.g., Becker et al. 2004; Morales-Rueda et al. 2006; 
Bramich et al. 2008), unveiling a variety of transients 
ranging from the common (Rau et al. 2008) to the unex- 
plained (e.g., Barbary et al. 2008). Untargeted ("blind") 
synoptic wide- field imaging surveys, such as the Texas 
Supernova Search (TSS; Quimby 2006) conducted with 
the ROTSE-III 0.45-m telescope (Akerlof et al. 2003), 
have uncovered a large number of rare transients, includ- 
ing the four most luminous supernovae (SNe) observed 
to date: SN 2005ap (Quimby et al. 2007), SN 2008am 
(Yuan et al. 2008a), SN 2006gy (Ofek et al. 2007; Smith 
etal. 2007, 2008b), and SN 2006tf (Smith etal. 2008a). 
Observations of these very luminous events are starting 
to allow the detailed physical study of the extrema in 
core-collapse SNe. They appear to be powered in part 
by their interaction with a highly dense circumstellar 
medium (CSM; see Smith et al. 2008a, and references 
therein), though other possibilities have been advanced. 
Clearly, the discovery of more such events would allow 
an exploration of the variety of the phenomenology as 
related to the diversity of progenitors and CSM. 

Recently, the TSS discovered yet another luminous 
transient on 2008 Apr. 26.23 (UT dates are used 
throughout this paper), which they suggested was a vari- 
able active galactic nucleus at a redshift z = 1.02 (Yuan 
et al. 2008b). Gezari & Halpern (2008) then hypoth- 
esized that the transient was a flare from the tidal dis- 
ruption of a star by a supermassive black hole. Miller 
et al. (2008) first identified SN 2008es as potentially 
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an extremely luminous Type II SN (see also Gezari 
et al. 2008a), and we later definitively confirmed this 
with further spectroscopic observations (Chornock et al. 
2008a); the event was assigned the name SN 2008es by 
the lAU (Chornock et al. 2008b). It is located at a = 
ll''56™49.06", 5 = -h54°27'24.77" (J2000.0). 

Here we present our analysis of SN 2008es, which 
is classified as a Type II-Linear (ILL) SN based on 
the observed linear (in mag) decline in the photomet- 
ric light curve (Barbon, Ciatti, & Rosino 1979; Doggett 
& Branch 1985). At z = 0.213, SN 2008es has a peak 
optical magnitude of My — —22.3, among SNe second 
only to SN 2005ap. Aside from the extreme luminos- 
ity, SN 2008es is of great interest since detailed ultra- 
violet (UV) through infrared (IR) observations provide 
a unique opportunity to study the mass-loss properties 
of an evolved post-main sequence massive star via its 
interaction with the surrounding dense CSM. A similar 
analysis of SN 2008es has been presented by Gezari et 
al. (2008b). 

The outline of this paper is as follows. We present 
our observations in §2, and the photometric and spectro- 
scopic analyses of this and pubHc (NASA) data in §3 and 
§4, respectively. A discussion is given in §5, and our con- 
clusions are summarized in §6. Throughout this paper 
we adopt a concordance cosmology of i7o = 70 km s~^ 
Mpc-\ = 0.3, and f^A = 0.7. 

2. OBSERVATIONS 

Here we present our ground-based optical and near- 
infrared (NIR) photometry and optical spectroscopy, 
along with space-based Swift UV, optical, and X-ray ob- 
servations. NIR photometry of SN 2008es was obtained 
simultaneously in J, H, and Ks with the Peters Auto- 
mated Infrared Imaging Telescope (PAIRITEL; Bloom 
et al. 2006) beginning 2008 May 16. To improve the 
photometric signal-to-noisc ratio (S/N), we stacked im- 
ages over multiple nights. For the Kg images the S/N 
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TABLE 1 

PAIRITEL Observations of SN 2008es 



TABLE 2 
KAIT Observations of SN 2008es 





Obs. window'' 


Filter 


Exp. time 




(day) 


(day) 




(sec) 




4602.24 


2.06 


J 


2895.91 


17.74±0.03 


4605.28 


6.09 


J 


3869.06 


17.68±0.03 


/I (^1 O 17 
40 iZ. i t 




J 


4606.78 




4615.15 


6.08 


J 


6458.90 


17.68±0.02 


4618.19 


2.19 


J 


6129.29 


17.77±0.03 


4632.17 


6.04 


J 


4944.24 


17.75±0.04 


4602.24 


2.06 


H 


2880.22 


17.62±0.06 


4605.27 


6.09 


H 


3751.34 


17.43±0.06 


4612.17 


2.09 


H 


4598.93 


17.59±0.06 


4615.15 


6.08 


H 


6451.06 


17.40±0.04 


4618.19 


2.19 


H 


6074.35 


17.53±0.05 


4632.17 


6.04 


H 


4865.76 


17.89±0.11 


Note. 


— PAIRITEL 


observations were 


stacked over 



multiple epochs to increase the S/N. 

*Mid-point between the first and last exposures in a single 
stacked image, reported as JD— 2,450,000. 

''Time between the first and last exposures in a single 
stacked image. 

''Observed value; not corrected for Galactic extinction. 

of the SN remained low, despite the stacks made over 
multiple epochs, and therefore we do not include these 
data in our subsequent analysis. Aperture photometry, 
using a custom pipeline, was used to measure the J and H 
photometry of the isolated SN calibrated to the 2MASS 
catalog (see Bloom et al. 2008). The resulting light 
curves are presented in Figure 1. The final PAIRITEL 
photometry is reported in Table 1. 

Optical photometry of SN 2008es was obtained in 
BVRI with the Katzman Automatic Imaging Telescope 
(KAIT; Filippenko et al. 2001) and the 1-m Nickel 
telescope at Lick Observatory beginning 2008 May 30. 
Point-spread function (PSF)-fitting photometry was per- 
formed on the SN and several comparison stars using the 
IRAF/DAOPHOT package (Stetson 1987) and transformed 
into the Johnson-Cousins system. Calibrations for the 
field were obtained with the Nickel telescope on three 
photometric nights. The final photometry from KAIT 
and the Nickel telescope is given in Tables 2 and 3, re- 
spectively. 

Additional optical photometry was obtained in 
UBVRI on 2008 June 02 and BVRI on 2008 Aug. 05 
with PFCam on the 3-m Shane telescope at Lick Ob- 
servatory. The data were reduced using standard tech- 
niques and aperture photometry was used to extract the 
SN flux. BVRI calibrations were done using the same 
stars as those used for the Nickel and KAIT observations. 
Calibrations for the U band were not obtained with the 
Nickel telescope. Therefore, we convert the Sloan Digi- 
tal Sky Survey (SDSS) colors of stars in the field to the 
U band with the color transformations of Jester et al. 
(2005), and use these stars for our iJ-band calibration. 
The final PFCam photometry is reported in Table 4. 

The Swift satellite observed SN 2008es during 13 
epochs between 2008 May 14 and Aug. 21. We 
downloaded the Ultraviolet/Optical Telescope (UVOT; 
Roming et al. 2005) data from the Swift data archive 
and analyzed the Level 2 sky image data in U, B, and V 
according to the photometry calibration and recipe by Li 
et al. (2006). The Swift UV filters (UVWl, UVMl, and 
UVW2) were reduced following Poole et al. (2008). The 
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17.67±0.05 


4624.72 




300.00 


17.81±0.05 


4626.72 




300.00 


17.88±0.04 


4628.72 




300.00 


17.88±0.05 


4630.72 
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17.90±0.07 


4632.74 




300.00 


18.01±0.06 


4634.77 




300.00 


17.93±0.07 


4636.73 




300.00 


18.01±0.08 


4639.70 




300.00 


18.24±0.11 


4643.74 




300.00 


18.20±0.09 


^Exposure 


mid-point, 


reported as 



JD-2,450,000. 
''Observed value; not corrected for Galac- 
tic extinction. 

final Swift UVOT photometry is reported in Table 5. 

Simultaneous observations of SN 2008es occurred with 
the Swift X-ray Telescope (XRT; Burrows et al. 2005), 
for which we confirm a nondetection of X-ray emission 
(see also Gezari & Halpern 2008). To place a limit on 
the source flux we assume a power-law spectrum with 
a photon index F = 2, absorbed by elements associated 
with Galactic H I. We took an extraction region of ra- 
dius 64 pixels (^2.5') and fit the PSF model around the 
ccntroid of the optical emission. We obtain a 3(t limiting 
flux of 9 X lO-i'"^ erg cm^^ g-i ^he 0.3 10 keV band 
for a total exposure of 54.7 ks. This represents an upper 
limit to the X-ray luminosity of SN 2008es of ~ 1.2 x 10^^ 
ergs s~^. 

We obtained spectra of SN 2008es on 2008 May 16.3, 
2008 May 29.3, and 2008 July 7.3 using the Kast spectro- 
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TABLE 3 
Nickel Observations of SN 2008es 
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Fig. 1.— Observed UV-optical-NIR light curves of SN 2008es. 
The data have not been corrected for Galactic or host-galaxy ex- 
tinction. We include ROTSE-IIIb unfiltered observations from the 
literature (open symbols; Gezari et al. 2008b), as wrell as our 
optical-NIR observations (filled symbols; KAIT, Nickel, PFCam, 
and PAIRITEL) and space-based UVOT observations from Swift 
(filled triangles). We adopt the discovery date of SN 2008es, 2008 
Apr. 26.23, as "day 0" for this SN. Low-order polynomial fits to 
each band have been overplotted to help guide the eye. 

graph on the Lick 3-m telescope (Miller & Stone 1993). 
Additional spectra were obtained on 2008 June 7.4 and 
2008 Aug. 3.3 using the Low Resolution Imaging Spec- 
trometer on the Keck I 10-m telescope (Oke et al. 1995) 
and on 2008 June 21.2 and 2008 June 23.2 using the R. 
C. Spectrograph on the Kitt Peak 4-m telescope, follow- 
ing the approval of our request for Kitt Peak Director's 
Discretionary Time. The spectra were extracted and cal- 
ibrated following standard procedures (e.g., Matheson et 
al. 2000). Clouds were present during several of the 
observations, making the absolute flux scales unreliable. 
The spectrograph slit was placed at the parallactic an- 
gle, so the relative spectral shapes should be accurate 
(Filippenko 1982), with the exception of the Kitt Peak 
spectra, which have a small amount of second-order light 
contamination at wavelengths longward of ^8000 A. The 
full SN 2008es spectral sequence is plotted in Figure 2. 



''Exposure mid-point, reported as 
JD-2,450,000. 

''Observed value; not corrected for Galac- 
tic extinction. 



TABLE 4 
PFCam Observations of SN 2008es 



^obs 


Filter 


Exp. time 


Mag'' 


(day) 




(sec) 




4619.81 


U 


1230.00 


17.73±0.07 


4619.83 


B 


630.00 


18.16±0.01 


4683.70 


B 


360.00 


21.39±0.15 


4619.83 


V 


930.00 


17.86±0.01 


4683.71 


V 


360.00 


20.19±0.05 


4619.79 


R 


930.00 


17.77±0.01 


4683.68 


R 


600.00 


19.67±0.04 


4619.85 


I 


1230.00 


17.63±0.01 


4683.70 


I 


360.00 


19.28±0.03 



''Exposure mid-point, reported as JD- 
2450000. 

Observed value; not corrected for Galac- 
tic extinction. 

The two Kitt Peak spectra show little evolution in the 
two days that separate them and have been combined 
to increase the S/N (day 33 in Figure 2). A log of our 
spectroscopic observations is presented in Table 6. 

We searched our spectra for the possible presence of 
narrow lines and were unable to positively identify any, 
either in emission or absorption. Therefore, without a 
detection of the host galaxy, we determine the redshift 
of SN 2008es directly from the SN spectrum. As the SN 



4 



Miller et al. 



TABLE 5 
UVOT Observations of SN 2008es 
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''Exposure mid-point, reported as 
JD-2,450,000. 

'^Observed value; not corrected for 
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Fig. 2.— Spectral evolution of SN 2008es. Spectra of SN 2008es 
(in color) are labeled with their ages in the rest frame of the su- 
pernova (z = 0.213) relative to the observed V-band maximum on 
May 13.3. The spectra become progressively redder as the SN ages 
and broad P-Cygni spectral features become more prominent with 
time. By +68 d, a broad emission feature of Ha is clearly present. 
The top spectrum (in black) is the earliest spectrum of SN 2005ap, 
the most luminous observed SN (Quimby et al. 2007). 

aged, broad P-Cygni spectral features appeared, includ- 
ing an emission feature near 7900 A that we identify as 
Ha near a redshift of 0.2. Other spectral features are 
consistent with a Type II SN at about that redshift. 

To get a more accurate redshift, we identified two sim- 
ilar reference spectra of the Type II-L SNe 1979C and 
1980K from the literature (Branch et al. 1981; Uomoto 
& Kirshner 1986). These spectral comparisons are shown 
in the bottom panel of Figure 3. We used the Super- 
Nova IDentification code of Blondin & Tonry (2007) to 
cross-correlate the day 68 spectrum of SN 2008es with 
the two reference spectra and derived a weighted-average 
z = 0.213±0.006, which we adopt throughout this paper. 
This redshift agrees with a fit to the broad Ha emission 
line in the day 68 spectrum, which yields a flux centroid 
of z = 0.210. The SN 2008es redshift determined using 
this method is correlated with the expansion velocity, but 
we assume that the bias due to this effect is small given 
that all three SNe have Ha emission lines of a similar 
width at late times. 

3. PHOTOMETRIC RESULTS 

Our dataset provides excellent broad-band coverage of 
SN 2008es from the UV to the NIR, which allows us to 
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TABLE 6 
Log of Spectroscopic Observations 



Age'' 


UT Date 


Instrument'^ 


Range 


Exp. Time 


Seeing 


Airmass 


Photometric? 


(days) 






(A) 


(s) 


(") 




(y/n) 
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2008-05-16.345 


Kast 


3300-8830 


1500 


2.4 


1.3 
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1.1 
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68 


2008-08-03.254 
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0.9 
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n 



''Age in rest-frame days relative to the observed V-band maximum on 2008 May 13.3. 

''Kast = Kast spectrograph on Lick 3-m telescope. LRIS = Low Resolution Imaging Spectrometer on Keck-I 10-m telescope. 
RC — R. C. Spectrograph on Kitt Peak 4-m telescope. 
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Fig. 3. — Spectral comparisons of SN 2008es to two SNe II-L. All 
spectra are labeled with their respective ages in rest- frame days rel- 
ative to maximum light, and prominent spectral features are iden- 
tified. The zero point of the flux scale is accurate for the SN 2008es 
spectra (in blue), while the comparison spectra (in black) are ver- 
tically offset. The top panel shows a comparison at an early epoch 
of SN 2008es to SN 1979C from 1979 Apr. 28 (Branch et al. 1981). 
The bottom panel shows a comparison at a later epoch to SN 1979C 
from 1979 May 28 (Branch et al. 1981) and SN 1980K (spectra 
from 1980 Nov. 15 and 17 combined; Uomoto & Kirshner 1986). 
Telluric absorption bands in the SNe 1979C and 1980K spectra are 
marked with a ® symbol. 

model changes in the spectral energy distribution (SED). 
In order to sample each photometric band onto a single 
set of common epochs, we fit lovif-order polynomials to 
each light curve, which we then interpolate onto a com- 
mon grid. NIR observations were only included on or 
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Fig. 4. — Photospheric and bolometric luminosity evolution of 
SN 2008es. Top: temperature evolution of SN 2008es based on BB 
fits (open circles) and inferred radius (open triangles). The error 
bars on the radius are likely to be slightly overestimated because 
the errors on the BB temperature and luminosity are correlated. 
Bottom: bolometric luminosity of SN 2008es derived via two in- 
dependent methods, BB modeling (closed squares) and integration 
of the total UV+optical (+NIR where available) flux (open dia- 
monds). The two methods agree to within < 20%. 

around epochs where we detected the SN. 

We create an SED at each of the common epochs and 
fit a single-component blackbody (BB) to the data fol- 
lowing the procedure described by Modjaz et al. (2008). 
Prior to the SED/BB fits we add a systematic term to 
the uncertainty in the photometric measurement in each 
band. This term is added in quadrature to the statistical 
uncertainty, and for the NIR (JH) is equal to 2%, in the 
optical (BVRI) we adopt 3%, while for U band we adopt 
10%, and the adopted UV {UVWl, UVM2, UVW2) sys- 
tematic uncertainty is 20%. Across all epochs of the BB 
fits the total per degree of freedom is 45.6/65, which 
suggests that our systematic terms may be slightly over- 
estimated. We assume no host-galaxy extinction (for fur- 
ther details see §5), and correct our measurements for the 
modest amount of Galactic reddening E{B — V) = 0.011 
mag (Schlegel et al. 1998). From the SED/BB fits we 
derive the temperature, radius, and luminosity of the SN 
as a function of time, as shown in Figure 4. The temper- 
ature and luminosity decrease while the radius increases 
with time, as expected for an expanding and cooling SN 
photosphere. 
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Fig. 5. — SEDs from SN 2008es corrected for Galactic extinction. 
We show two representative SEDs of SN 2008es at '^27 and ~65 
rest-frame days after discovery. The SEDs are well fit by a single 
BB component. Note that relative to the model BB there is excess 
flux in the UVW2 band in both epochs, while the UVM2 band 
shows excess in the later epoch. All UV measurements that exhibit 
a clear excess relative to a BB are shown in gray and are excluded 
from the fits. 

In order to determine the time of V-band maximum 
light, which at the redshift of this SN roughly corre- 
sponds to rest-frame B, we convert the early time g, r, 
and i' photometry taken with the Palomar 1.5-m tele- 
scope from Gezari et al. (2008b) to VRI using the 
color equations from Jester et al. (2005). Following a 
quadratic fit to these data, we find that the observed 
y-band maximum occurred on ~2008 May 13, '^15 rest- 
frame days after discovery. The observed peak is V — 
17.8 mag, which corresponds to My — —22.3 mag, with 
a scatter about our fit of 0.04 mag. 

In Figure 5 we show two representative fits to the SED 
of SN 2008es at 26.6 and 65.4 rest-frame days after dis- 
covery. In all epochs we observe excess emission relative 
to a single BB in the bluest of the Swift filters, UVW2, 
while starting around 55 rest-frame days after discovery 
there is excess emission in both UVW2 and UVM2. This 
excess was also observed in the Type Il-Plateau (II-P) 
SN 2006bp, where it was attributed to complex line blan- 
keting by Fe-peak elements (Immler et al. 2007). The 
blue excess is readily identified in the UV color curves 
(e.g., UVW2-UVWI) of SN 2008es, which evolve to- 
ward the red until ^^50 rest-frame days after discovery, 
at which point the UVW2~UVWl and UVM2-UVWI 
colors become progressively more and more blue. These 
points were excluded from our SED/BB fits. To confirm 
that we were not underestimating the bolomctric lumi- 
nosity of the SN, we directly integrated the flux in the 
SED using the same method as Modjaz et al. (2008). 
We show the results of this direct integration in the 
lower panel of Figure 4, and in all cases we find that 
the BB model and direct integration method agree to 
within <20%. 

The observed linear decline in the light curve of 
SN 2008es leads us to classify it as a Type II-L SN. 
We measure a bolomctric photometric decay rate of 
0.042 mag d~^, which is slightly faster than the F-band 
(roughly rest-frame B) rate of 0.036 mag d~^. This rate 
is also slightly slower than the rest-frame i?-band de- 
cay of both SN 1979C (0.046 mag d-^ Panagia et al. 



1980) and SN 1980K (0.055 mag d"^ Barbon, Ciatti, & 
Rosino 1982). Furthermore, integrating the bolomctric 
light curve from 15 to 83 rest-frame days after discovery 
yields a total radiated energy of ~9 xlO^" ergs, compa- 
rable to the canonical 10^^ ergs deposited into the kinetic 
energy of a SN. If we assume the same bolomctric correc- 
tion factor to the observed V'-band light curve both pre- 
and post-maximum, and we include the early data from 
Gezari et al. (2008b), we find the total radiated energy 
of SN 2008es over the first 83 d after discovery is ^^1.1 
xlO'^^ ergs. 

The three Swift UVOT observations taken more than 
80 d after the discovery of SN 2008es show tentative ev- 
idence for a significant reduction in the SN photomet- 
ric decline rate (see Figure 1). We note, however, that 
the late-time UVOT measurements have significant er- 
rors, and these observations may be consistent with the 
early trend seen in the light curve. Gezari et al. (2008b) 
suggested that these observations were evidence that the 
SN luminosity made a transition to being powered by ra- 
dioactive heating from ^^Co. With only two observations 
of this decline, we caution against prematurely identify- 
ing the ^^Co decay tail, and instead argue that the ob- 
servations are inconclusive at this time. For instance, 
it would be possible to mimic the behavior of ^^Co de- 
cay with a late-time light echo, or interaction with an 
extended CSM. We note that SN 2006gy also showed 
evidence for a "flattening" several months after explo- 
sion (Smith et al. 2007), but further observations taken 
the following year indicated that this was not clearly the 
^^Co tail (Smith et al. 2008b). Our last optical spectrum 
(Fig. 2) also shows no evidence that the ejecta were be- 
coming optically thin on that date, as would be expected 
if SN 2008es were transitioning to the radioactive decay 
tail. Late-time photometry will be necessary to conclu- 
sively identify if and when the light curve of SN 2008es 
made a transition to being powered by ^^Co decay, and 
thus allow a measurement of the amount of ^^Ni synthe- 
sized in the explosion. 

4. SPECTROSCOPY 

The SN 2008es spectral sequence plotted in Figure 2 
is labeled with ages (in the rest frame) relative to the 
observed V-band maximum in order to facilitate com- 
parison with SNe 1979C and 1980K in Figure 3. Our 
first two spectra of SN 2008es (at -f3 and -1-13 d rela- 
tive to maximum light) show a smooth and featureless 
blue continuum with no identifiable spectral features. In 
particular, we do not detect an emission feature near 
5650 A (4660 A in the rest frame) seen in earlier spectra 
of SN 2008es taken between 2008 May 1 and 2008 May 
8 (Yuan et al. 2008b; Gezari et al. 2008b). Gezari et 
al. (2008b) attribute this feature solely to He II A4686, 
but we note that transient emission features seen in the 
early spectra of some SNe Iln at similar wavelengths are 
due to a combination of He II A4686 and the Wolf-Rayet 
C III/N III A4640 blend (Niemela, Ruiz, & Phillips 1985; 
Leonard et al. 2000). 

In addition, the spectrum of the extremely luminous 
SN 2005ap from — 3 d is shown at the top of Figure 2. 
The SN 2005ap spectrum is also very blue, but has addi- 
tional spectral features that are not seen in our SN 2008es 
spectrum from a similar epoch. Quimby et al. (2007) 
identified the strongest spectral feature in SN 2005ap, 
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the "W" -shaped feature near 4200 A, as being due to a 
blend of C III, N III, and O III with an expansion velocity 
of about 20,000 km s^^. 

Our next spectrum, at +21 d, is noticeably redder 
and is the first to show strong spectral features. Both 
SN 1979C at a similar epoch and SN 2008es show a blue 
continuum with weak, low-contrast lines of H and He I 
lines present mostly in absorption (Fig. 3). However, Ha 
is present only in emission and is weak in SN 2G08es. The 
Ha line in SN 2008es has a full width at half-maximum 
intensity (FWHM) of 10,000 km s^^ and an equivalent 
width of only 22 A. (Both of these values have large un- 
certainties due to difficulties in defining the continuum 
for a line with such a low amplitude.) One difference 
between the two objects is the lower apparent velocities 
in SN 2008es. The H/3 and He I A5876 lines of SN 1979C 
have absorption minima at velocities of 9700 and 8900 
km s-^ In SN 2008es, these values are 6000 and 5700 
km s~^, respectively, although we caution that the ex- 
act values are correlated with the assumed value for the 
redshift. 

Over the next two months, the spectra of SN 2008es 
plotted in Figure 2 became redder, reflecting the cool- 
ing photospheric temperature evolution discussed above. 
In addition, the spectral features first seen in the day 
-|-21 spectrum gradually became more prominent. The 
SN features are still muted in amplitude relative to those 
expected in a normal Type II SN. This may be an exam- 
ple of the "top-lighting" effect described by Branch et al. 
(2000), where continuum emission from interaction with 
CSM illuminates the SN ejecta from above and results in 
a rescaling of the amplitudes of spectral features. 

By the time of our day -1-68 spectrum, the P-Cygni 
spectral features due to H Balmer lines, Na I, and 
Fe II become more prominent and the overall appearance 
starts to resemble that of normal SNe II. The Ha profile 
lacks an absorption component, which may be common 
to SNe II-L and not SNe H-P (e.g., Schlegel 1996, Fil- 
ippenko 1997). The broad Ha emission extends (at zero 
intensity) from -9000 to 9000 km s-\ with a FWHM 
of 9500 km s^^. This velocity width is intermediate be- 
tween that of the SNe 1979C and 1980K spectra (FWHM 
w 10,600 and w 8200 km s-\ respectively). 

Unlike in most core-collapse SNe, the velocity of the 
minimum of the H/3 line increased over time, from 6000 
km s^^ at day 21, to 8700 km s~^ at day 68, as shown 
in the bottom panel of Figure 6. The exact values of the 
velocity depend directly on the assumed redshift, but 
the trend is independent of those uncertainties. Unfor- 
tunately, the other absorption lines are mostly blended 
(except for Ha, which does not show an absorption 
component), so we cannot isolate the velocity trend in 
other spectral features without spectral modeling. The 
top panel of Figure 6 shows the evolution of the He I 
A5876/Na I A5892 blend. At early times He I dominates 
the blend, and as the ejecta cool at later epochs Na I 
should dominate, resulting in an ~800 km s~^ redward 
shift of the rest wavelength. After taking into account 
that redward shift, the top panel of Figure 6 also shows 
some weak evidence for a blueward shift of the absorption 
minimum to ~8000 km s~^. 

The only supernova known to us to show increasing ab- 
sorption velocities over time is the peculiar SN lb 2005bf. 




Velocity (1000 km s"') 

Fig. 6. — Velocity evolution of absorption minima. The spectra 
are labeled with dates after maximum light in the same manner as 
in Figure 2. In the bottom panel the apparent blueshift of the H/3 
absorption minimum increases over time. A vertical dashed line 
at —6000 km s~^ marks the velocity of the absorption minimum 
on day 21 to guide the eye. In the top panel, the evolution of the 
He I A5876/Na I A5892 blend is plotted, with A5876 used as the 
zero point of the velocity scale. The vertical dashed line at —5700 
km s^^ marks the absorption minimum on day 21 to guide the eye. 

The trend of increasing absorption velocities was only 
visible in the three He I lines, but not in other lines such 
as Ca II H&K (Modjaz 2007). Tominaga et al. (2005) 
explained the effect as being due to progressive outward 
excitation of He I by radioactive ^^Ni in the interior as 
the ejecta expanded and the density decreased, an effect 
that seems to be of little relevance to SN 2008es. A more 
likely possibility is that blending with some unidentified 
line is shifting the velocity of the apparent absorption 
minimum. At late times H/3 dominates its region of the 
spectrum, but at early times He II A4686 could possibly 
be contributing emission to the blue wing of the H/5 ab- 
sorption profile. Another scenario is that some unusual, 
but as yet unidentified, radiative transfer effect is affect- 
ing the wavelength of the apparent absorption minimum 
in SN 2008es. For example, if the optical depth in an op- 
tically thick line increases sufficiently rapidly with time 
(e.g., as the ejecta cool and recombine), then the appar- 
ent absorption minimum could move blueward (Jeffery 
& Branch 1990). None of these suggestions explain why 
such an effect would be present only in SN 2008es and 
not in normal SNe. 

5. DISCUSSION 
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Fig. 7. — Rest-framo brightness evolution of five of the most 
luminous known SNe. For SN 2008es wc derive the absolute mag- 
nitude using UVOT and the KAIT/Nickel V-band magnitudes 
(filled green squares) and the KAIT/Nickel i?-band magnitudes 
(red squares), assuming that day is the discovery date. Early 
g-band photometry from Gezari et al. (2008b), converted to V- 
band using the color equations of Jester et al. (2005), is also shown 
(open green squares). The light curve for SN 2005ap (orange) is 
derived from unfiltered photometry in Quimby et al. (2007), and 
SN 2006gy (gray circles) is also unfiltered from Smith et al. (2007). 
ij-band photometry of SN 2006tf (blue circles) is from Smith et al. 
(2008a), and r/r' photometry of SN 2005gj (dotted line) is from 
Prieto et al. (2007). The light curve for SN 2006tf is shifted by 
+16 d from that in Smith et al. (2008a); since the explosion date 
is not known, we chose to align its time of peak luminosity with 
those of SNe 2005ap and 2008es. 

5.1. The Physical Nature of SN 2008es 

The photometric evolution of SN 2008es (see Fig. 7) 
is much faster than that of other very luminous SNe like 
SNe 2006gy (Smith et al. 2007; Ofek et al. 2007), 2006tf 
(Smith et al. 2008a), and 2005gj (Prieto et al. 2007), 
and its spectrum also betrays no evidence for the strong 
CSM interaction seen in these other SNe Iln, in the 
form of narrow lines from the CSM or intermediate- width 
Ha from the post-shock shell. Indeed, considering both 
its photometric and spectroscopic evolution, we suggest 
that SN 2008es is most like the overluminous SN 2005ap 
(Quimby et al. 2007) and thereby most closely resembles 
a Type II-L that is 4-5 mag more luminous than typical 
SNe II-L (Richardson et al. 2002). 

To power the tremendous luminosity of SN 2008es 
with radioactive decay would require an initial ^^Ni 
mass of --10 M© (see, e.g., Smith et al. 2007). This 
^^Ni mass would very likely need to be generated in a 
pair-instability explosion (Barkat, Rakavy, & Sack 1967; 
Bond, Arnett & Carr 1984), but this ^^Ni mass seems 
problematic given that it would be larger than the mod- 
est envelope mass indicated by the relatively fast rise 
and decay time (see below). In addition, the photomet- 
ric decline of 0.042 mag d~^ is faster than that of '^^Co, 
0.0098 mag d^^, making radioactive heating unlikely as 
the dominant source of the energy. 



Despite the lack of a Type Iln spectrum, the most 
likely interpretation seems to be that the high luminos- 
ity of SN 2008es is the result of converting shock energy 
into visual light. This can be accomplished, in principle, 
if the shock kinetic energy is thermalized throughout a 
massive envelope, like a normal SN II-P or II-L, but with 
a much larger initial radius of (2-3) xlO^^ cm (based on 
the peak luminosity and the evolution in Fig. 4). Were 
the initial radius much smaller than this it would prove 
difficult to convert ^10^^ ergs of kinetic energy to '^10^^ 
ergs of radiation because there would be significant adi- 
abatic losses. The apparent temperature evolution from 
~15,000 to 6000 K over a time period of ~66 d (Fig. 4), 
reminiscent of other normal SNe II, suggests that the 
recombination photosphere is receding through a cool- 
ing envelope in SN 2008es. In this scenario, the usual 
narrow/intermediate- width Ha emission that is taken as 
a signpost of CSM interaction might be avoided if the 
CSM shell is initially very opaque, and if the shock en- 
counters no further CSM at larger radii (see Smith et al. 
2008a). 

Since the required initial radius exceeds that of the 
largest known red supergiants (see Smith et al. 2001) 
by a factor of 20-30, it requires the envelope to be an 
unbound, opaque CSM shell ejected prior to the SN ex- 
plosion instead of a traditional bound stellar envelope. 
Similar models were suggested for SN 2005ap (Quimby 
et al. 2007), SN 2006gy (Smith & McCray 2007), and 
SN 2006tf (Smith et al. 2008a), implying CSM envelope 
masses of 0.6, ^10, and 18 Mq, respectively. The corre- 
sponding CSM mass for SN 2008es would be roughly '^5 
Mq in this scenario, because its evolution and expansion 
speeds are slower than those of SN 2005ap. This CSM 
mass of ~5 Mq allows the observed Ha Doppler veloc- 
ities to remain faster than in SNe 2006gy and 2006tf, 
where the heavy CSM shells decelerated the shocks to 
only 4000 and 2000 km s^^, respectively (Smith et al. 
2007; Smith et al. 2008a). The smaller CSM mass for 
SN 2008es relative to SNe 2006gy and 2006tf would also 
explain the comparatively short rise and decay observed 
for SN 2008es, as the light diffusion time for this SN is 
much shorter than in SNe 2006gy and 2006tf (Smith & 
McCray 2007). The putative envelope ejection preceding 
SN 2008es must have occurred 10-100 yr prior to the ex- 
plosion (for an unknown progenitor wind speed of wcsm 
= 10-100 km s~^), indicating a progenitor mass-loss rate 
of order 0.05-0.5 Mq yr~^. This is much larger than any 
steady stellar wind (see Smith et al. 2007, and refer- 
ences therein) , providing another case of impulsive mass 
ejection in the decades immediately preceding some SNe. 

An alternative analysis of SN 2008es, submitted for 
publication shortly after this paper was made available 
electronically, also suggests that the extreme luminos- 
ity is powered via interaction with CSM (Gezari et al. 
2008b). However, the initial version of Gezari et al. 
(2008b) argues against an unbound CSM shell ejection 
prior to the SN explosion, and instead argues that the 
CSM is created by a dense progenitor wind. 

5.2. The Host of SN 2008es 

Currently there is no conclusive detection of the host 
galaxy of SN 2008es, so its metallicity and correspond- 
ing implications for the pre-SN evolution are not known. 
From SDSS DR6 images (Adelman-McCarthy et al. 
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2008) , wc derive a Sa upper limit of m^' > 22.7 mag 
at the SN position, which translates to roughly My > 
—17 A mag at the SN redshift (neglecting any K correc- 
tions). Thus, the putative host galaxy is significantly 
less luminous than an galaxy and could be compara- 
ble to or fainter than the Small Magellanic Cloud (with 
My = -16.9 mag)-*. 

Without a definitive detection of the host galaxy we 
assume that there is no host extinction. This assump- 
tion is further supported by both the low luminosity of 
the host galaxy and our BB fits (see §3) that show ex- 
cellent agreement with our UV measurements, where the 
effects of host extinction would be especially prominent. 
In analogy with the underluminous hosts of SNe 2006tf 
(Smith et al. 2008a) and 2005ap (Quimby et al. 2007), 
deep imaging after the SN has faded may uncover the 
host galaxy nearly coincident with the SN position, which 
would allow further constraints to be placed on extinction 
in the host. These three extreme SNe and their underlu- 
minous hosts may be hinting that very luminous SNe II 
preferentially occur in low-luminosity host galaxies. 

5.3. Rates of Extremely Luminous SN 2008es-like 

Events 

Over the past three years the TSS has successfully 
found tens of SNe with a surprising rate of unusual 

objects. The list of such SNe includes SN 2005ap, 
SN 2006gy SN 2006tf, SN 2008am, and now SN 2008es. 
However, this apparent high anomaly rate is probably the 
result of combining a huge survey volume with an intrin- 
sically rare class of objects. Following Quimby (2008), 
we compare the volume probed by the TSS for SNe la, 
and bright core-collapse SNe. The highest redshift of the 
~30 SNe la found by the TSS is z f« 0.1, while SN 2005ap 
was found at 2 « 0.3. The comoving volume scanned by 
the TSS is therefore ~23 times bigger for SN 2008es-like 
objects than for SNe la. Of the five bright objects found 
by the TSS, only two, SN 2005ap and SN 2008es, can 
be considered together as a class. Thus, the comparative 
rate is about 30/2 x 23 ~ 350 times smaller. 

The KAIT SN search (FiUppenko et al. 2001) has dis- 
covered about 400 SNe la in the past 10 years in tar- 
geted nearby galaxies, hence the expected number of 
SN 2008es-like objects is of order 1. This is consistent 
with the nondetection of any such SN. 

Using the light curve of SN 2008es and the accumu- 
lated observations of about 10, 000 spiral galaxies in the 
KAIT sample which were observed regularly during the 
past ten years, we find an upper limit on the rate (per 
unit ii'-band luminosity) of such SNe to be about 1/160 
times the local SN II rate (Leaman, Li, & Filippenko 

2009) . Furthermore, we note that if SN 2008es-hke events 
preferentially occur in small, low-metallicity galaxies (see 
§ 5.2), the KAIT sample, which targets large, nearby 
galaxies, would be biased against the discovery of such 
a SN. While there are several SN search surveys with 
larger search volumes than the TSS, such as the SDSS- 
II Supernova Survey (Sako et al. 2008), the full details 
of these surveys (cadence, total survey volume, candi- 
date SN color cuts, etc.) are not available at this time, 

A galaxy ~9" to the NE of the SN position is unlikely to be 
the host itself even if found to be at a similar redshift as SN 2008es; 
this would require significant massive star formation at a projected 
physical distance of ~31 kpc. 



and therefore we cannot directly compare them to the 
TSS. We conclude that the number of luminous Type II- 
L SNe, as discovered by the TSS, does not appear to be 
the result of a statistical fluke, and predict many such 
detections with future wide-field synoptic surveys. 

6. CONCLUSIONS 

We have reported on our early observations of 
SN 2008es, which at a peak optical magnitude of My = 
—22.3 is the second most luminous SN ever observed. We 
argue that the extreme luminosity of this SN was likely 
powered via strong interaction with a dense CSM, and 
that the steep decline in the light curve, 0.042 mag d~^, 
indicates that the radioactive decay of ^^Co is likely not 
the dominant source of energy for this SN. Integration of 
the bolometric light curve of SN 2008es yields a total ra- 
diated energy output of > 10^^ ergs. The optical spectra 
of SN 2008es resemble those of the luminous SN 1979C, 
but with an unexplained increase in the velocity of the 
II/3 absorption minimum over time. 

We also examined the rate of discovery of extremely 
luminous SNe by the Texas Supernova Search and find 
that their discovery of the five most luminous observed 
SNe in the past four years is probably not a fluke; several 
more such detections are expected in the coming years. 

Finally, what behavior can we expect from SN 2008es 
at late times, roughly 1 yr or more after discovery? Re- 
gardless of whether the peak luminosity was powered by 
radioactive decay or optically thick CSM interaction (see 
Smith et al. 2008b), a SN can be powered by strong 
CSM interaction at late times if the progenitor had a 
sufficiently high mass-loss rate in the centuries before 
exploding. We have seen examples of both: SN 2006tf 
had strong Hq emission indicative of ongoing CSM in- 
teraction at late times (Smith et al. 2008a), whereas 
SN 2006gy did not (Smith et al. 2008b). SN 1979C was 
less luminous at peak than those two SNe, but it has 
been studied for three decades because its late-time CSM 
interaction is powering ongoing emission in the radio, op- 
tical, and X-rays (Weiler et al. 1981; Fesen & Matonick 
1993; Immler, Pietsch, & Aschenbach 1998). With such 
an extraordinarily high peak luminosity, a late-time IR 
echo such as that seen in SN 2006gy (Smith et al. 2008b) 
is also likely if SN 2008es has dust waiting at a radius 
of ~0.3 pc. Alternatively, if SN 2008es were powered in 
whole or in part by radioactivity, a large mass of ^^Ni 
should be evident in the late-time decline rate. 
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